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The Ising spin-glasses are investigated on three dual pairs of hierarchical lattices, using exact 
renormalization-group transformation of the quenched bond probability distribution. The goal is 
to investigate a recent conjecture which relates, on such pairs of dual lattices, the locations of the 
multicritical points, which occur on the Nishimori symmetry line. Towards this end we precisely 
determine the global phase diagrams for these six hierarchical spin-glasses, using up to 2.5 x 10 9 
probability bins to represent the quenched distribution subjected to an exact renormalization-group 
transformation. We find in all three cases that the conjecture is realized to a very good approxima- 
tion, even when the mutually dual models belong to different spatial dimensionalities d and have 
different phase diagram topologies at the multicritical points of the conjecture and even though the 
contributions to the conjecture from each lattice of the dual pair are strongly asymmetric. In all 
six phase diagrams, we find reentrance near the multicritical point. In the models with d — 2 or 
1.5, the spin-glass phase does not occur and the phase boundary between the ferromagnetic and 
paramagnetic phases is second order with a strong violation of universality. 
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I. INTRODUCTION 



II. THE CONJECTURE 
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The phase diagram structure of spin-glasses remains an 
open field of inquiry, since most approaches to the prob- 
lem rely on approximations. Any exact analytical result 
in this area is thus very valuable, both for the direct infor- 
mation it provides and as a test for approximation meth- 
ods. Over the last few years striking progress has been 
made combining the replica method, duality, and sym- 
metry arguments Q, 0, M 0, , an approach which has 
yielded the exact locations of the multicritical points in 
the Ising and Potts spin-glasses on the square lattice and 
in the four-dimensional random-plaquette gauge model. 
The most recent result in this series |5j is a general con- 
jecture relating the multicritical point locations of any 
spin-glasses on a pair of mutually dual lattices. In sup- 
port of the conjecture, estimates based on Monte Carlo 
simulations were given for Ising spin-glasses, in d = 2, 
on the dual pairs of triangular and hexagonal lattices 
and, in d = 3, on the dual pairs of bilinear and lattice- 
gauge interactions on the cubic lattice. In both cases, 
within the numerical limitations, the conjecture is ap- 
proximately satisfied. 

We propose here to extensively test the conjecture in 
an alternative fashion using hierarchical lattices [|| 0, , 
by loo king at Ising spin-glasses on mutually dual pairs 
0, El El 0| of such lattices. These constitute ideal 
testing grounds, since an exact renormalization-group 
transformation for the quenched bond probability dis- 
tribution can be constructed for such lattices, yielding 
global phase diagrams and critical properties. Accord- 
ingly, the location of the phase boundaries and of the 
multicritical points are precisely determined. We thus 
investigate three pairs of hierarchical lattices, and in the 
end find that the conjecture is very nearly satisfied for 
all of them. 



The Ising spin-glass is given by the Hamiltonian 



-f3H = Jas,s 

(ij) 



3 > 



(1) 



where Sj = ±1 at each site i, (ij) denotes a sum over 
nearest-neighbor pairs of sites, and the bond strengths 
Jij are equal to + J with probability 1 — p and — J with 
probability p. The limits p — and p — 1 correspond to 
purely ferromagnetic and purely antiferromagnetic sys- 
tems respectively. 

To give a brief overview of the conjecture, let us con- 
sider the model on an arbitrary lattice, and treat the 
randomness through the replica method, where the sys- 
tem is replicated n times and the n — > limit is even- 
tually taken, in order to get results for the physical sys- 
tem. The partition function of the n-replicated system 
after averaging over randomness, Z n , can be expressed 
entirely as a function of n + 1 "averaged" combinations 
of edge Boltzmann factors, e ±J , associated with nearest- 
neighbor bonds 0.FJ|. These averaged Boltzmann factors, 
Xk(p, J), k = 0, . . . ,n, have the form 

x k (p, J) = pe- (n - k)J e kJ + (1 - p)e {n ~ k *> J e~ kJ , (2) 

where the ktYi factor corresponds to a configuration with 
a parallel-spin bond in n — k replicas and an antiparallel- 
spin bond in k replicas 0. Thus, 



Z n (x (p,J), xi(p, J), x n (p, J)) 



(3) 



The partition function on the dual lattice, Z*, can be 
expressed in a similar form, 



K = Z ni x o(Pi J),Xi(p,J),..., x* n {p, J)) 



(4) 
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with the dual counterparts to the averaged Boltzmann 
factors given by 



_ t j\ n— 2k / _ j j \ 2k 

e + e \ ( e — e 



V2 



V2 

ri-2fc-l 



(5) 



_j Jx 2k+l 



V2 

for < 2k < 2k + 1 < n. Z n and Z* are related as 

Z n (x (p, J), . . . , x n (p, J)) 

= 2 a Z:(x*(p,J), ...,<(p,J)), (6) 

where a is a constant, which can be eliminated by us- 
ing Eq. 10 evaluated at two different sets of parameters, 
(pi, Ji) and (p2, J2), giving a relationship of the form 



Z n (x (pi, Ji), . . . , x n (pi, Ji)) 

' Z*( x o(P2, J2), Xn(P2, J2)) 

' Z n {x* Q (p2, J2), Xn(P2, J2)) ■ 



(7) 



The individual partition functions Z n can be rewrit- 
ten by extracting xq, the averaged Boltzmann factor 
corresponding to an all-parallel spin state, thus effec- 
tively measuring the energy of the system relative to this 
state Q: 

Z n (x ,xi, ...,x n ) = x„ B Z n {u 1 ,u 2 , ...,u n ), (8) 

where Nb is the number of bonds in the lattice, and the 
reduced variables are w, = Xi/xo- Eq. Q becomes 

[x (pi, Ji)x Q (p 2 , J 2 )] NB Z n (u 1 (pi, Ji), . . . , u„(pi, Ji)) 
■ Z*(u 1 (p 2 ,J 2 ), • • ■ , u n (p 2 , J 2 )) 

- ^l)^0 (P2, J 2 )] ArB 2*(ui(pi, Ji), . . . , <(pi, Ji)) 

' ^n(Wl(P2, J2), ■ ■ ■ , <(P2, J2)) • 

(9) 

In general, the form of Eq. © is too complicated to 
yield useful information relating the locations of phase 
transitions. However, the multicritical points in both 
original and dual systems are expected to lie 0, Q, 
on the Nishimori line []], which simplifies the relation. 
Furthermore, the conjecture advanced in Ref. [f| states 
that, for the multicritical points (pi m , Jim) of the origi- 
nal system and (p2m, J2m) of its dual, Eq. (|5J) is satisfied 
when the leading Boltzmann factors xq from each side 
are equal, 

Xo(Plm, Jlm)xo(P2m, -hm) = X* {pi m , Jl m )xl(j>2m, -hm) ■ 

(10) 

Since (j>i m , Jim) and (p2m, he on the Nishimori line, 



1 -Plr 
Pirn 



= 2J 2 „ 



1 ~ P2r, 
P2m 



From Eqs. (0 and JSJ , Eq. ljlfj|l gives 

(P?r + (1-Plm) n+I )(p^ 1 + (1-P2 m ) n+1 ) = 2-. (12) 

Finally taking the limit, n — > 0, one obtains the condition 



H{pim) + H(p 2 



(13) 



where H(p) = —plog 2 p — (1 — p)log 2 (l — p). As ex- 
pressed in Eq. {15). the conjecture is asserted to hold for 
multicritical points of Ising spin-glasses on any pair of 
mutually dual lattices p|. 



III. THE MULTITUDE OF ISING 
SPIN-GLASSES ON HIERARCHICAL LATTICES 



(11) 



Hierarchical lattices |6|,l7|,|8j are constructed by replac- 
ing every single bond, in a connected cluster of bonds, 
with the connected cluster of bonds itself, and repeating 
this step an infinite number of times. These provide mod- 
els exactly solvable by renormalization group, with which 
complex problems have been studied and understood. 
For example, frustrated [lfij , sp in-glass [l7j . random- 
bond and random-field |l9j| . Schrodinger equation 
pof , lattice- vibration [2l| , dynamic sc aling [23] , aperiodic 
magn et [23j . complex phase diagram |24j . and directed- 
path |25| systems, etc., have been solved on hierarchical 
lattices. 

To test the conjecture of Eq. l(T5|) . we study Ising spin- 
glasses on the dual pairs of hierarchical lattices, depicted 
in Figs.^El and|5J Each lattice in a given pair is the dual 
of the other. These particular choice of lattices are mo- 
tivated by their properties under renormalization-group 
transformation as related to physical lattices. The hierar- 
chical lattices of Fig. 1(a) and (b) yield the two variants 
of the Migdal-Kadanoff recursion relations ptil |2^ | for 
dimension d — 2 with length rescaling factor 6 = 3. Sim- 
ilarly, the lattice in Fig. 2(a) yields a Migdal-Kadanoff 
recursion relation for d = 3, 6 — 3. Its dual lattice in 
Fig. 2(b) has d = 1.5,6 = 9. (The two variants of the 
Migdal-Kadanoff recursion relations correspond to mu- 
tually dual hierarchical lattices only in d = 2.) Lastly, 
the hybrid lattice in Fig. 3(a) is interesting because it has 
been shown to give very accurate results for the critical 
temperatures of the d — 3 isotropic and anisotropic Ising 
model JiU]. This lattice has d = 3, 6 = 3, while its dual 
in Fig. 3(b) has d = 1.5,6 = 9. 



IV. EXACT RENORMALIZATION-GROUP 
TRANSFORMATION OF HIERARCHICAL 
SPIN-GLASSES 

For a pure system, the renormalization-group trans- 
formation on a hierarchical lattice consists of a decima- 
tion by summing over the internal sites in each of the 
connected clusters making up the lattice (the right-hand 
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(b) 




FIG. 2: Another pair of mutually dual hierarchical lattices. 
The Migdal-Kadanoff recursion relations are exact for lattice 
(a) with d — 3, 6 = 3. Its dual lattice, in (b), has d = 1.5, 
6 = 9. 



sides of Figs. 11130 . Thus, the hierarchical lattice con- 
struction process is reversed, as each connected cluster is 
replaced by a single renormalized bond. The decimation 
can be expressed as a mapping, 



J i>j> 



(14) 



where the set {Jij} are all the bonds within the connected 
cluster of the original system and J[,j, is the renormalized 
bond between sites i 1 and j' of the rescaled system. In 
the pure case, all Jy bonds are independent for ij, and 
the implementation of Eq. I|14[l is straightforward. 

When quenched randomness is added to the system, 
the renormalization-group transformation is expressed in 
terms of quenched probability distributions [18j, where 
the quenched probability distribution V\J[, j,) in the 
rescaled system is calculated from P(Jij) in the original 
system through the convolution 



* J 



S (4,-, - R({Jij})) 



Here the product runs over all the bonds ij in the con- 
nected cluster of the original system between sites i' and 



The recursion of the quenched probability distribu- 
tion, Eq. (| 1 51) . is implemented numerically. The prob- 
ability distribution is represented by histograms, each 
histogram being specified by a bond strength and an as- 
sociated probability. Thus, for the spin- glass problem, 
the starting distribution consists of two histograms, one 
at J with probability 1 — p, and one at — J with probabil- 
ity p. Eq. Ijl5|l dictates the convolution of 9 probability 
distributions for the lattices of Fig. ^ an d the convolu- 
tion of 27 distributions for the lattices of Figs.[21and|3| In 
this task, computational storage limits can be maximally 
exploited by factorizing Eq. Ijl5(l into an equivalent series 
of pairwise convolutions, each of which involves only two 
distributions convoluted using an appropriate R function. 
The types of pairwise convolutions needed are a "bond- 
moving" convolution, with 



Rbm{Ji 1 j 1 , Jiij-i) ~ Jiljl "b Ji2j2 

and a decimation convolution, with 



1 



hi 



(16) 



(17) 



which is just the standard decimation transformation for 
a two-bond Ising segment. 

Consider the hierarchical lattice in Fig. da). If Pmit 
is the initial probability distribution, a series of pair- 
wise convolutions which yields the total convolution of 
Eq. (|15fl for this lattice is: (i) a bond-moving convolution 
of 7-mit with itself, yielding V\ ; (ii) a bond- moving con- 
volution of V\ with Vyoxt , yielding Vi ; (iii) a decimation 
convolution of Vi with itself, yielding V3 ; (iv) a decima- 
tion convolution of Vz with Vi-, yielding T'Wai- For the 
lattice in Fig. Ufa), the series is: (i) a decimation convo- 
lution of Pinit with itself, yielding V\ ; (ii) a decimation 
convolution of V\ with Vi ni t, yielding Vi\ (hi) a bond- 
moving convolution of V2 with itself, yielding V3; (iv) a 
bond-moving convolution of V3 with itself, yielding V4; 
(v) a bond-moving convolution of V4 with itself, yielding 
V5] (vi) a bond-moving convolution of V5 with V%, yield- 
ing "Pfinai. For the lattice in Fig.^a), the series is: (i) a 
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bond-moving convolution of V- m it with itself, yielding V\ ; 
(ii) a decimation convolution of V\ with itself, yielding 
Vi\ (hi) a decimation convolution of V% with Vi, yield- 
ing V3] (iv) a bond- moving convolution of V3 with itself, 
yielding V4] (v) a bond-moving convolution of V4 with 
itself, yielding V5 ; (vi) a decimation convolution of Vmit 
with itself, yielding Ve ; (vii) a decimation convolution of 
Ve with Pinjt, yielding 7-V; (viii) a bond- moving convolu- 
tion of 7-V with V5, yielding PsixaX- As for the dual lattices 
in Figs. [Jb), |21[b), and|3fb), the series of pairwise convo- 
lutions are identical to their counterparts above, except 
that each bond-moving is replaced by a decimation, and 
vice versa. 

Since the number of histograms that constitute the 
probability distribution increases rapidly with each 
renormalization iteration, a binning procedure is used 
when the desired (large, namely up to 2.5 x 10 9 ) num- 
ber of histograms is reached: Before every pairwise con- 
volution, the histograms are placed on a grid, and all 
histograms falling into the same grid cell are combined 
into a single histogram in such a way that the average and 
the standard deviation of the probability distribution are 
preserved. Histograms falling outside the grid, represent- 
ing a negligible part of the total probability, are similarly 
combined into a single histogram. Any histogram within 
a small neighborhood of a cell boundary is proportion- 
ately shared between the adjacent cells. In the current 
study, the binning procedure is done separately for J > 
and J < 0. After the convolution, the original number of 
histograms is reattained. 

In the current study, 40,000 bins are generally 
used, representing the renormalization-group flows of 
80,000 variables, requiring the calculation of 40,000 
local renormalization-group transformations at each 
renormalization-group iteration. The numerical results 
converge rapidly with increasing bin number. For 
maximal accuracy in determining the exact locations 
of the multicritical points, in the immediate vicin- 
ity of these points we used at least 1,000,000 his- 
tograms, representing the renormalization-group flows 
of 2,000,000 variables, requiring the calculation of 
1,000,000 local renormalization-group transformations at 
each renormalization-group iteration. It should thus be 
noted that our analysis is an exact numerical solution of 
Ising spin-glasses on hierarchical lattices. 



V. RESULTS 

Global phase diagrams for the various hierarchical lat- 
tices are obtained from the renormalization-group flows 
of the probability distributions. Each phase has a cor- 
responding sink, namely a completely stable fixed distri- 
bution. The boundaries between phases flow to unstable 
fixed distributions, analysis of which yields the order of 
the phase transition and the values of the critical expo- 
nents of second- and higher-order transitions. All the 
phase diagrams are plotted in terms antiferromagnetic 




0.1 0.2 0.3 0.4 0. 

Antiferromagnetic bond concentration p 

FIG. 4: Phase diagrams for the two hierarchical lattices in 
Fig. 1, with the solid lines indicating second-order phase tran- 
sitions between the ferromagnetic (F) and paramagnetic (P) 
phases. In each diagram the multicritical point, separating 
two different types of second-order boundary, is marked by a 
dot, and the Nishimori symmetry line is drawn dashed. The 
phase diagrams were calculated with 40,000 probability bins, 
except for the vicinity of the multicritical points, where for 
higher precision 1,000,000 probability bins were used. 



bond concentration p versus temperature 1/J. The dia- 
grams are symmetric around p = 1/2, with the ferromag- 
netic phase in the p < 1/2 half-space mapping onto the 
antiferromagnetic phase in the p > 1/2 half-space. Thus 
in the figures only the p < 1/2 portions are shown. 

Fig. 4(a) and (b) show the phase diagrams for the dual 
pair of hierarchical lattices in Fig. 1(a) and (b) respec- 
tively. The phase structure of both diagrams is topolog- 
ically identical to that of the d = 2 Ising spin-glass on 
a square lattice, which is only natural considering that 
the d = 2, b = 3 Migdal-Kadanoff recursion relations are 
exact on these hierarchical lattices. The p = transi- 
tion temperatures of the two models are related by the 
duality algebra [2^ 



sinh(2 J lc ) sinh(2 J 2c 



(18) 



which is also true for the two other pairs of mutually dual 
hierarchical models. Furthermore, the p = transition 
temperatures in Fig. 4(a) and (b) are related 

US by 



-1 T-l 

J 2c ' 



(19) 
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0.1 0.2 0.3 0.4 0. 

Antiferromagnetic bond concentration p 




0.1 0.2 0.3 0.4 0. 

Antiferromagnetic bond concentration p 



FIG. 5: Phase diagrams for the two hierarchical lattices in 
Fig. 2, with the solid lines indicating second-order phase 
transitions between the ferromagnetic (F), paramagnetic (P), 
and spin-glass (SG) phases. In each diagram the multicrit- 
ical point is marked by a dot, and the Nishimori symmetry 
line is drawn dashed. The phase diagrams were calculated 
with 40,000 probability bins (250,000 bins for the inset in the 
top figure), except for the vicinity of the multicritical points, 
where for higher precision 10 6 and 2.5 x 10 9 probability bins 
were used in (a) and (b) respectively. 



FIG. 6: Phase diagrams for the two hierarchical lattices in 
Fig. 3, with the solid lines indicating second-order phase 
transitions between the ferromagnetic (F), paramagnetic (P), 
and spin-glass (SG) phases. In each diagram the multicrit- 
ical point is marked by a dot, and the Nishimori symmetry 
line is drawn dashed. The phase diagrams were calculated 
with 40,000 probability bins (250,000 bins for the inset in the 
top figure), except for the vicinity of the multicritical points, 
where for higher precision 2.5 x 10 7 and 2.25 x 10 8 probability 
bins were used in (a) and (b) respectively. 



since the mappings of the interaction constant in the 
repetition of renormalization-group transformations dif- 
fers only by an initial bond strengthening by a factor of 
6 d ~ 1 ; note that Eq. i fHty does not apply to < p < 1, 
since there the bond-moving is not a mere multiplica- 
tive strengthening, but a (b )-fold convolution of the 
probability distributions that alters this distribution in 
a non-simple way. Eq. 1|19[) is also not applicable to the 
two other pairs of mutually dual models, since the rep- 
etition of renormalization-group transformations are not 
differentiated by only a preliminary bond-moving. 

In each of Fig. 4(a) and (b), a ferromagnetic phase at 
low temperatures and low p is separated from the dis- 
ordered paramagnetic phase by two second-order phase 
boundaries, meeting at a multicritical point. (In a nar- 
row neighborhood of all multicritical points in our re- 
sults, reentrance is observed: paramagnetic, then ferro- 
magnetic, then paramagnetic or spin-glass phases are en- 
countered as temperature is lowered at fixed p.) The two 
second-order boundaries flow to distinct unstable proba- 
bility distributions with different critical exponents, con- 



stituting a strong violation of universality |17| and con- 
sistent with the prediction, generally, of the absence 
of first-order transition under quenched randomness in 
d = 2. [3(j As expected from symmetry considerations, 
the multicritical points fall 0, Q, 0] precisely on the 
Nishimori line [| as seen in Table I. As also seen in Table 
I, H{pi m ) + H{p2m) = 1.0172, so that the conjecture is 
realized to a very good approximation. 

Fig. 5 shows the phase diagrams for the dual pair of 
hierarchical lattices in Fig. 2. While Fig. 5(b) has the 
same phase topology as the diagrams in Fig. 4, being at 
d = 1.5 below the spin-glass lower-critical dimension, a 
different structure occurs in Fig. 5(a). Here the d = 3, 
b = 3 Migdal-Kadanoff relations are exact on the hierar- 
chical lattice, and for low temperatures in the vicinity of 
p = 1/2 there exists a spin-glass phase. The multicriti- 
cal point occurs where the ferromagnetic, paramagnetic, 
and spin-glass phases meet. As expected both multi- 
critical points lie directly on the Nishimori line. From 
Table I we see that H(p lm ) + H(p 2m ) = 0.9829, so that 
the conjecture is realized to a very good approximation, 
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Figure 


Plm, Jim 


J- 1 


P2m, J 2m 


j- 1 


H(p lm ) 


H(p 2m ) 


H(jpim) + H(p 2m ) 


Fig. 4 


0.1735,1.2810 


1.2810 


0.06620, 0.7557 


0.7557 


0.6656 


0.3516 


1.0172 


Fig. 5 


0.1851, 1.3494 


1.3494 


0.05128, 0.68546 


0.68545 


0.6911 


0.2918 


0.9829 


Fig. 6 


0.2473, 1.7973 


1.7973 


0.02796, 0.5636 


0.5636 


0.8070 


0.1840 


0.9911 



TABLE I: Locations of the multicritical points in the phase diagrams of Figs. 4-6 (corresponding to the hierarchical lattices 
of Figs. 1-3) . is the value calculated from p% m using Eq. (f 1 1 pi for the Nishimori line and turns out equal to , for both 
i = 1 and 2. The quantities H(pi m ) that enter the conjecture and their sums are also given. 



even when the mutually dual models belong to differ- 
ent dimensionalities d and have different phase diagram 
topologies at the multicritical points of the conjecture. 

The phase diagram structures in Fig. 6, correspond- 
ing to the dual pair of hierarchical lattices in Fig. 3, are 
similar to those of Fig. 5, illustrating dimensions above 
and below the spin-glass lower-critical dimension. Again 
the multicritical points for both cases lie directly on the 
Nishimori line. In this case H(pi m ) + H(p2 m ) = 0.9911, 
and the conjecture is realized to a very good approxima- 
tion, again for mutually dual models belonging to differ- 
ent dimensionalities d and having different phase diagram 
topologies at the multicritical points of the conjecture. 

Thus, we find that for all three mutually dual pairs of 
hierarchical lattices, the conjecture relating the locations 
of the multicritical points is satisfied to a very good ap- 
proximation. This is all the more remarkable, since, as 
seen in Table I, the contributions of H(pi m ) and H(p2m) 
to the conjecture are strongly asymmetric. However, it 
should be noted that (1.0172,0.9829,0.9911), while being 
very close to 1, arc different from integer 1. In our numer- 



ical implementation of the convolutions of the probability 
distributions, the results have converged to the precision 
of the digits shown in Table I. Further increase of the 
already very large number of probability bins does not 
change the entries in the table. Further tests of the con- 
jecture, using other systems, would be very useful. Simi- 
lar to our current study, the use of hierarchical lattices to 
study phenomena linked to mutually dual lattices, e.g., 
Ref. [31|, would also be very useful. 
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